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Abstract 

We calculate the two loop long range effect on the proton decay effective 
Lagrangian. Numerical calculation for suppression factor gives ^4i(2-loop) = 
0.321 for the value of the strong coupling constant a s (mz) — 0.116. Two loop 
effect to more general effective Lagrangian is also given. 

1 Introduction 

Proton decay is one of the most important prediction of grand unified theory 
(GUT). The main decay mode p — > K + + V in the minimal supersymmetric (SUSY) 
SU(5) GUT has been searched for with the underground experiments, and the lifetime 
bounds of this mode have been given by Kamiokande and 1MB as r > 1.0 X 10 32 
years and r > 0.62 x 10 32 years, respectively. [l] These bounds are rather close to the 
theoretical prediction [2] [3] [4] in case of the SUSY breaking parameters of less than 1 
TeV and the Superkamiokande is expected to improve these bounds drastically in the 
near future. It is now desirable to make a more careful estimation of the theoretical 
bounds. 

In the following we estimate the two loop QCD correction to the effective proton 
decay Lagrangian. This correction is estimated by solving the renormalization group 
equation from the weak energy scale to the low energy scale, say 1 GeV. 
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In order to formulate the renormalization group equation we need to estimate the 
two loop renormalization constant of the effective Lagrangian, which is not calculated 
so far though the two loop effect of the quark wavefunction, the quark masses and 
(3- function are calculated in the literature [5]; this is why reliable calculation with the 
two loop approximation is not made yet. In this paper, the renormalization constant 
is calculated for rather general effective Lagrangian for proton decay in the two loop 
approximation and the renormalization group equation of the coefficient of the proton 
decay effective operator is given. The numerical estimation of the effect is also made. 

2 The renormalization constant of the effective 
operator 

The dominant proton decay effective Lagrangian at the weak scale is written as [6] 

£ eff = C^e^iq^qfMJi) + h.c, (1) 

where q and £ are left-handed Weyl spinors of quark and lepton field respectively. 
Here a,/?, 7 are color indices, i 7 j 7 kj are flavor indices and the bracket ( , ) means 
the inner product for Weyl spinors. 

We shall also consider the proton decay Lagrangian with mixed quark chiralities, 

like 

£' eff = DWea^tcfi, o^m^i) + h - c -> (2) 

or 

£" eff = E^ kl e a(h (¥t^MJi) + h.c, (3) 

to discuss the two loop corrections to different four fermion vertices. Here a M denotes 
the Pauli matrices and q c is the Weyl spinor for antiquark field. 

The long range effect of renormalization on the coefficient C^ kl in Eq.(l) is to be 
evaluated by the ratio of its value at lGeV to that at the weak scale mz — 90GeV, 
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that is, 
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In literatures, however, the following quantity is more frequently used: 
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where \i represents the renormalization point and m u or means the up-quark or 
the down-quark mass respectively. 

The one loop correction has already been given as [2] 
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Here the strong coupling constant at the renormalization point \i — mz{mz — 90GeV) 
is taken as a s (mz) — 0.116. 

We take account of only QCD correction since the strong coupling constant is 
much larger than the weak or electromagnetic coupling constants. Up to the two loop 
level, there are six different types of Feynman diagrams for the QCD correction as 
given in the Fig.'s 1(a), 1(b), • • •, 1(f). 

Each type of diagram in Fig.l is calculated by straightforward calcutations and 
the result is given in Table 1. In the table, MS scheme is adopted [7], 

^ = ^-7 + hi(47r), (7) 

where n — 4 — 2e ( n-dimensional space-time ) and 7 is an Euler number. 

In Table 1, we take A = gll^nf = a s /(47r) ( g s is a strong coupling constant ) 
and 77 as the gauge parameter in the gluon propagator in the following form: 
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From Table 1, we obtain the renormalization constant Zq as 

A A 2 ( 74 8N F 29 3 9 

A 2 ( 2N F 31 7 2 \ 

including the one loop effect. Here A = for the left-handed ( or right-handed ) 
quarks and lepton as in Eq.(l), while A = —10/3 for the mixed fermion chiralities as 
in Eq.(2), Eq.(3). 

3 The renormalization group equation of C^ kl 

The quark wavefunction renormalization constant Z qi QCD /^-function /?, quark 
mass renormalization factor 7 m are calculated in the literature[5] as followes; 

A 9 A 2 9 { (a 2 25\ 3 9 ) 

Z q = l--aT 2 --T 2 \[ J+ a + -}c 2 -N F t--T 2 ^ 

A 2 9 { ( a 2 3a\ a 2 9 ] . . 

+ w T 2 {(t + tJ C2 + tH' (10) 

P{9s) = -^(yC 2 -^)-^(^C 2 2 -4 W T 2 -fwC 2 ),(ll) 

lm {9s) = 8A + 8A 2 {^- ^j. (12) 

with N F being the number of quark flavors, a — 1 — 77, C 2 — 3, T 2 = | and t — |. 
The running of g s and the quark mass m q are given by 
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Vrr m q = -lm{g s )m q . (13) 



Combining Eq.(9) and Eq.(lO), we get the renormalization group equation of the 
coefficient C njkl as 

a/i I a/i J 



= - |<L4 + A 2 (y + ^p + A) + 0(A 3 )} C ijkl . (14) 

In the following we take A = unless otherwise stated, since the dominant decay 
mode in the minimal SUSY SU(5) GUT is derived from the Lagrangian in Eq.(l) with 
left-handed quarks and lepton. 

Solving Eq.(14), one can get the two loop long range effect on the proton decay 
effective Lagrangian between weak scale and lGeV: 
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The result of numerical calculation is given in Table 2. Two loop correction changes 
Al, but Ri does not change very much. Uncertainty of the strong coupling constant 
has a great influence on these corrections. Here we give some comments on the 
Lagrangians with mixed chiralities. We find the one loop correction is common to all 
types of the local four fermion proton decay Lagrangian with three quarks and one 
lepton, while we find that the two loop corrections to them are common except for 
the diagram of the type Fig. 1(e) as seen in the Table 1. 

4 Conclusion 

We have obtained the two loop correction for the long range effect on the proton 
decay effective Lagrangian and got the result i?i(2-loop) = 1.33 and ^(2-loop) = 
0.32 for a s (mz) — 0.116. We have found that the two loop correction of Ai is 
appreciably large though the two loop four fermion vertex correction itself has a 
small contribution compared with the quark mass and gauge coupling renormalization 



correction. We note that the four fermion vertex correction, though found to be 
small, is necessary to obtain a consistent evaluation of the two loop renormalization 
effect. We have also found that the two loop correction to the effective proton decay 
Lagrangian depends on the combination of the chiralities of the four fermion fields 
while the one loop correction does not. 
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Figure caption 

Figure 1 : The types of diagrams in the two loop level. The bold line in the center 
represents that we take the inner product of and qj and that of qk and 1%. Fig. 1(a) 
and Fig. 1(b) represent a typical two loop diagram with a quark and a gluon self- 
energy correction subgraph, respectively. Fig. 1(c) is that with a gluon-quark vertex 
correction subgraph. Fig. 1(d) contains a four fermion vertex correction subgraph. 
Fig. 1(e) has gluon crossing propagators. In Fig. 1(f) gluons are emitted from three 
quark lines and meet at the same point. 

Table captions 

Table 1 : Contribution from each type of diagrams in Fig.l to Zq- Definitions of 
A{— a s /(47r)), e, r\ ( gauge parameter ), Np and A are given in the text. ( A = for 
the left-handed ( or right-handed ) quarks and lepton as in Eq.(l), while A = —10/3 
for the mixed fermion chiralities as in Eq.(2), Eq.(3). ) 

Table 2 : Numerical results for Ri and A^. The value of the strong coupling 
constant is a s {mz) — 0.116, 0.111 and 0.121 for (a), (b) and (c), respectively. 
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Figure 1: The types of diagram in the two loop level. The bold line in the center 
represents that we take the inner product of qi and qj and that of q^ and l\. Fig. 1(a) 
and Fig. 1(b) represent a typical two loop diagram with a quark and a gluon self- 
energy correction subgraph, respectively. Fig. 1(c) is that with a gluon-quark vertex 
correction subgraph. Fig. 1(d) contains a four fermion vertex correction subgraph. 
Fig. 1(e) has gluon crossing propagators. In Fig. 1(f) gluons are emitted from three 
quark lines and meet at the same point. 
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type 


number of diagrams 


coefficient of A 2 je 2 


coefficient of A 2 je 


a 


6 


- - V )(2 - 77) 


16 , 

-j (1-1) 


b 


3 


2N F 3 

5 1 — 77 

3 2 ' 


26 8N F 7 3 2 

1 1 — 77 77 

3 9 2 ' 4 ' 
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6 


52 25 2 
y- 1777 + ^ 


112 2 
9 9^ 
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16 8 
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3 3 ' 
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1 








total 


28 


2N F 31 7 2 

25 ?7 + -?7 2 

3 2 ' 2 ' 


74 8N F 29 3 9 

h + —?? ?? + A 

3 9 9' 4 ' 



Table 1: Contribution from each type of diagrams in Fig.l to Zc- Definitions of 
A{— a s /(47r)), e, 77 ( gauge parameter ), N F and A are given in the text. ( A = for 
the left-handed ( or right-handed ) quarks and lepton as in Eq.(l), while A = —10/3 
for the mixed fermion chiralities as in Eq.(2), Eq.(3). ) 
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(a) a s {mz) — 0.116 
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A L 


1-loop 


1.31 


0.449 


2-loop 


1.34 


0.321 


(b) a s (m z ) = 0.111 
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2-loop 


1.31 
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(c) a s {m z ) = 0.121 






A L 


1-loop 


1.33 


0.421 


2-loop 


1.38 


0.277 



Table 2: Numerical results for Rl and Al- The value of the strong coupling constant 
is a s (mz) — 0.116, 0.111 and 0.121 for (a), (b) and (c), respectively. 
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